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ABSTRACT 

in ■ 

^\ ' For over 50 years attempts have been made to explain the properties of nuclear 

T-H , matter in terms of constituent nucleons with very little success. Here we will in- 

pH ' vestigate one class of many possible models, string-flip potential models, in which 

, flux-tubes are connected between quarks (in a gas/plasma) to give a minimal over- 

' all field configuration. A general overview of the current status of these models, 

, along with some of our recently finished work, shall be given. It shall be shown that 

' these models seem promising in that they do get most of the bulk properties of 

I , nuclear matter correct with the exception of nuclear binding. Finally we will con- 

^ ' elude with a brief discussion on ways to improve the string-flip potential models in 

, an attempt to obtain nuclear binding (currently we are investigating short range 

0^ ' one-gluon exchange effects - some preliminary results shall be mentioned). 

m ■ 

in- 
i/-^ , 1. Overview 

5 ■ 

r-| ■ The main objective of our work is to attempt to describe nuclear matter in terms 

P-i! of its constituent quarks. A difficult task indeed, for over the past 50 years many 

Q^l attempts have been made with very little success. The main difficulty is due to the 

nonperturbative nature of QCD. The most rigorous method for handling multiquark 
systems to date is lattice QCD, but given the magnitude of our problem it appears 
• ^ ! unlikely to be useful in the near future, due to its computationally intensive nature. 

^ I As a result we must consider more phenomenological means. 

' The basic idea here is to construct models which are motivated by lattice QCD 

theory and nucleon based models of nuclear matter. A very crude model should be 
able to get most of bulk properties of nuclear matter correct: i.e. 

• Nucleon gas at low densities with no van der Waals forces. 

• Nucleon binding at higher densities. 

• Nucleon swelling and saturation of nuclear forces with increasing density. 

• Quark-gluon plasma at extremely high densities. 

The are many models out there that attempt to fill this shopping list but we have 
found none that covers it completely. 
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Here, we shall nestrict ourselves to a particular class of models, called string-flip 
Potential Models,tlu'@ which fills all the items on our list with the exception of nucleon 
binding. The later part of our discussion shall be a way of perhaps improving these 
models in an attempt to get nucleon binding. 



2. String-Flip Potential Model 

Sting-flip potential models are models which are inspired by lattice QCD in that 
they attempt to mimic the fiux-tube dynamics. In the most general setting, one 
assumes that the quarks, in question, move slowly enough that their fields have enough 
time to reconfigure themselves in order to minimize the overall potential energy: i.e. 

V = mm{ J2 v(?m •••?„) I U {Qm ■ ■ ■ Qn} = {qi ■ ■ ■ Qn,}} , (1) 

{qm---qn} {m...n} 

where the A''^ quarks are placed in a cube of side L and subjected to periodic boundary 
conditions, to simulate continuous quark matter. The sum is over all gauge invariant 
sets {qm-'-Qn} of quarks, such that at least one element from each set lies inside 

a common box, whose disjoint union, U , makes up the complete colour singlet set 
{qi---qNq} of Nq quarks. It is easy to see that this potential allows for complete 
minimal quark clustering separability at low densities without suffering from van der 
Waals forces. 

At present these models are quite crude in that they do not include short range 
one-gluon exchange phenomena and spin effects, and are flavour degenerate. 

2.1. SUi{3) & SUh{3) Models 

To simplify matters we shall assume that the multiquark potential, eq. runs 
over triplets of quarks and that the colour is fixed to a given quark: i.e. 

V = mm{ ^ v(rr,fg,rb)| \J {qrqgqb} = {qi ■ ■ ■ qN,}} ■ (2) 

{qrqgqb} {rgb} 

For a linear model, SU£{3), the three body potential, v(fj., fg, Fb), is given byi 
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(3) 



where Fj^ = r, — r,- , ^^gj, = (r^g + rgjj + r^j.)/3 , and A is the area enclosed by the triangle 
Argb . 



For a harmonic oscillator model, SUh{3), our three body potential is simply 



Vh(rr,rg,fb) — |A;^rgb 




(4) 



where we have assumed quarks of equal mass (c/. ref. |i|). 

The linear model is inspired by lattice QCD whereas the harmonic oscillator is used 
because it has some nice analytical properties which allows us to make a consistency 
check. 

2.2. Computation 

Here we use a variational procedure to find the binding energy per nucleon, Eb 
(= E{p) — E{0)), and the correlation length, /3, as a function of nucleon density, p. 
The variational wave function, 



^a/3p — -^Correlation '^'S'/ater ; (5) 

is a trial function which consists of the product of a symmetric correlation piece. 



Correlation length 

I'rms , 
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E (/^^rgb) 
{rgb} 



Variational parameter 
1.75 ifSU£(3) 
2 if SU,,(3) 
and is found by 
minimizing E(/3 — 0). 



(6) 



over the set of triplets 
{rgb} which minimizes V. 



and an antisymmetric Fermi piece. 
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Product of Slater determinants containing 

= sin(^ Hi ■ ij + Si) , 
where 6i — or Tr/2 , and (ni)a(=a;,y,2) = 0, ±1, ±2, . . . 

This gives a smooth interpolation between a correlated (nucleon) gas at low densities 
and an uncorrelated (Fermi) gas at high densities. 

The results of the variational procedure are shown in Fig. 1. The solid lines show 
the minimal Eb{p) and /3(p) curves, where we have minimized E{p,j3). The integrals 
involved in evaluating E{p, (3) were done using the Metropolis algorithmiH along with 
a bag of tricks to minimize cpu time.§ 
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Fig. 1: EeiMeV), (a), and f3{fm-^), (b), verses p{fm-^) for 5?i(3) (★) and SUnii) (o), where 
mq = imMeV/c^, a = QMCeV/fm and k k. ?,2AAM eV / f rri^E T\ve dashed curves represent 
the remnants of the minimal Eb{p,0) and Eb{p,(3) trajectories. 

From Fig. l.a we see that there is an overall saturation of nucleon forces followed 
by a transition (the kink in the solid line) to a Fermi gas, but no nucleon binding. 
Fig. l.b illustrates an overall swelling of nucleons with increasing p until they become 
infinite beyond the transition point. 



3. The Show So Far! 

We have demontrated that string-flip potential models gets the overall bulk prop- 
erties of nuclear matter correct with the exception of binding. 

So what can we do to obtain binding? 

a) What about allowing the colour to move around? 

— This seems unlikely; it has been investigated in an SUh{2) model by HP.i 

b) What about allowing higher order flux-tube topologies? 

— This seems unlikely; it been investigated in an - albeit slightly ill|^ - SUh{S) 
chain model by HP.Q 

c) What about "a)" and "b)"? 

— This seems^ unlikely; it has been investigated by many six quark (non-string- 
fiip) modelsQO which suggests six quark states hke to dissociate into two nucle- 
ons. 

d) What about including short range forces? 

— This seems likely; we are currently investigating this possibility. 



3.1. Short Range Forces 

At the moment we are investigating the inclusion of colour Coulomb effects. The 
difficulty arises here when trying to combine perturbative and nonperturbative fields 
in a consistent fashion. To simplify our study we consider an SUi{2) potential model, 
which looks a lot like SU^i?)) {cf. ref. |), for qq pairs: i.e. 



0- r,; 



■ 1 



if Tij > To 



asXij if < ro 




(8) 



Coulomb 



-3/4, 1/4 for unlike and like colours respectively, and as = 0.1. Fig. 2 



where Xij ■ 
illustrates the dynamics of the model. 



Fig. 2: Consider configuration (a) of quarks, with r > ro, about 
to move to (b), s.t. two of them are within r < rg. Then the 
procedure is to draw a bubble of ro away from the two, (b), and 
to cut the flux-tubes at the boundary and insert virtual qq pairs, 
(c). Once the potential is computed the configuration is restored 
to (b) before the next move is made. 



Nonperturbative 
Piece 



a) 







Cut here 



\ Perturbabtive 
P I bubble 



Virtual qq pair insertion 



V ~ ar + X\ 



Notice that this model allows us to construct colourless objects because of the inser- 
tion the virtual qq pairs. These virtual quarks are used as a tool to get the overall 
length of fiux-tube correct. They are not used in computing the Coulomb term how- 
ever, as the field energy is already taken into account by the "real" quarks in the 
bubble. In general, once the bubbles have been determined, one must reconfigure the 
fiux-tubes in order to minimize the linear part of the potential. 





Fig. 3: Eb{p), (a), and f3{p), (b), for ro Ofm (- -) and ro = O.lfm (— ). 



Some preliminary results of these models are shown in Fig. 3. For ro = Ofm we get 
back the same results - as expected - and for ro = 0.1 fm we get a slightly different 
answer. 

Although the model is currently for SUe{2) it should be easy to extend it to a 
full SUe{3) model with all the one-gluon exchange phenomena. The results presented 
here are very preliminary as our wave function is slightly sick. By this we mean that 
we used our old wave function, eq. where the correlation piece was over the 
modified flux-tubes and a new piece was thrown in to take care of the local Coulomb 
interactions as they occurred: i.e. 

\T/ ^, ^XLincar) (Coulomb) ^(Fermi) /q\ 

However because the fiux-tubes and bubbles can now be created or destroyed the wave 
function is, in general, no longer smooth and continuous: i.e. we are not guaranteed 
a variational lower bound. Currently we are working to rectify the matter and hope 
to have some results out in the near future. 



4. Closing Remarks 

We have demonstrated that string-flip potential models appear to be fairly good 
candidates for explaining nuclear matter at the constituent quark level. Crude models 
seem to get most of the bulk properties with the exception of nuclear binding. We have 
also shown that it is fairly reasonable to assume that these models must be extended 
if we are ever going to achieve binding. In this light we have modified the model to 
include one-gluon exchange effects. Currently it only includes Coulomb effects, which 
we hope to have results for sometime soon. Further down the round we would like to 
include the rest of the one-gluon exchange phenomena. 
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